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O , 

O , Abstract 

Tj" ■ In loop quantum cosmology, a difference equation for the wave function describes 

the evolution of a universe model. This is different from the differential equations 
that arise in Wheeler-DeWitt quantizations, and some aspects of general properties 
of solutions can appear differently. Properties of particular interest are boundedness 
and the presence of small-scale oscillations. Continued fraction techniques are used 
', to show in different matter models the presence of special initial conditions leading 

to bounded solutions, and an explicit expression for these initial values is derived. 

> 

u ■ 1 Introduction 



Quantum gravity can in general be expected to provide a foundation which is radically 
different from classical smooth background geometries in that space and time should be 
discrete. Evolution equations are then no longer given by differential equations in space 
and time, but by difference equations which in general can be approximated by differential 
equations only on large scales. This is realized explicitly in loop quantum cosmology 
[T], based on loop quantum gravity El El, where one obtains difference equations for 
the evolution of a universe [Sj. This has physical implications, in particular concerning 
the singularity issue, but also conceptual ones since solutions of difference equations can 
generically behave differently from those of differential equations. 

In this paper we focus on properties of solutions of isotropic loop quantum cosmology 
in the large volume region, where one expects differential equations to provide a good 
approximation [HI U\. There are several cases where the boundedness of solutions is of 
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interest, for instance when a probability interpretation is used jH] or the physical inner 
product is studied jHj- From the continuum limit one can draw generic expectations, as in 
quantum mechanics, since one only needs to know whether regions of interest are classically 
forbidden. In such a case, one expects exponential behavior of solutions of a differential 
equation and bounded solutions only in exceptional cases. However, since these issues 
can depend sensitively on precise values of a solution, a direct analysis of the difference 
equation without detour through the continuum limit is preferable. 

This analysis is performed here with a general confirmation of the expectations. As a 
spin-off of the investigation we obtain a semi-analytical procedure which allows to deter- 
mine initial values for bounded solutions. The resulting formula is much stronger than a 
numerical trial and error procedure, and does not only apply to those few models where 
exact solutions are known. 

2 Isotropic loop quantum cosmology 

Dynamical properties of a universe in isotropic loop quantum cosmology fOJ El E] are 
determined by a wave function ip^ subject to a difference equation of the type 



in terms of volume eigenvalues = (7£||/z|/6) 3//2 ^3] which also show the geometrical 
meaning of the label A precise model is specified by the choice of a matter field 

and Hamiltonian ^ ma tter (/•*) and by fixing k to be zero for a flat model or one for a model 
of positive spatial curvature. The parameter 7 > is the Barbero-Immirzi parameter of 
loop quantum gravity > which fixes the scale of spatial discreteness. There are other 

versions of the equation, for instance obtained after a symmetrization procedure which 
changes the coefficients but not their large fi behavior. Our analysis will also apply to 
this equation. However, sometimes equations of higher order appear for which a similar 
analysis will be more complicated, as commented on later in this paper. Similarly, partial 
difference equations obtained for anisotropic models ^E] require more care. 

At small \i the difference equation can be used to understand quantum aspects of the 
classical singularity at \i = [13 UHl HE], while at large one expects semiclassical 
behavior Here, we are interested in asymptotic properties of solutions at large in 
particular their boundedness or oscillatory properties. For simplicity, we base the analysis 
on phenomenological matter systems which do not require an additional matter field. In 
those cases, the equation takes the asymptotic form 



with s M := e 1 ^/ 4 ?/^ and a(/i) 2 := 7£p|/i|/6. Matter Hamiltonians of interest are the cos- 
mological constant term H(a) = Aa 3 , dust with H(a) constant, radiation H(a) oc a -1 , or 




(1) 



s M+4 - (2 + kj 2 )^ + s^ A 



(2) 
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a free massless scalar with plane wave initial conditions where H(a) oc a~ 3 . In all these 
cases, the right hand side of (J2J) is of the form — A\fi\ v where v is constant in a given model, 
but can in general take both signs. We simply denote constants characterizing the matter 
content by A, even if they do not correspond to a cosmological constant. Furthermore, 
without loss of generality we drop the constant k and define 7 to be zero in flat models 
(after having absorbed it in A on the right hand side). This brings us to the general 
equation 



to be studied. 

Since the step size of the difference equation is four, defined for a wave function on the 
real line, values of ip^ between any /x and /1 + 4 are unrelated by the equation. Depending 
on the initial conditions, there can thus be oscillations on arbitrarily small scales, which 
for intuitive reasons one would like to be suppressed in semiclassical regimes. If this is the 
case, a solution is called pre-classical [IH1IZ|- This is partially a matter of initial conditions, 
but even within a given sequence at an integer lattice \x + 4Z can small scale oscillations 
develop. Whether or not this can be suppressed depends now on the difference equation 
and thus has to be decided by an asymptotic analysis. Oscillations on larger scales are 
allowed and show the presence of a classically allowed regime. Also the occurrence of these 
oscillations can be seen by a general analysis. 

In regimes where there are no large scale oscillations, which then are classically forbid- 
den, one expects exponential behavior of solutions. Generically, the exponentially growing 
contribution will dominate such that the issue of boundedness of solutions arises in partic- 
ular in models showing a classical recollapse [Sj. 

Both issues, pre-classicality and boundedness, have been used for an understanding 
of the semiclassical limit [Hj and the physical inner product |5], even though the precise 
status of the conditions is not clear. Nevertheless, in those investigations one often needs 
to find special solutions, such as bounded ones, which are difficult to detect numerically in 
a reliable manner. It is thus necessary to have analytical tools available. For oscillations, 
generating function techniques have been used [23 Ell - but this is applicable only in special 
cases (such as integer v in our equation, and with ignoring the absolute value signs). The 
methods used here are more general and also apply to the boundedness issue. 

It may happen that additional conditions will force /i to take discrete values [22] (enu- 
merated by but not necessarily identical to integer numbers). It is therefore reasonable 
to construct tools in such a way that this possible discretization will be built in. To do 
so, we first notice that any real number can be represented as \x = n + 5 with n G Z and 
< 5 < 1, and define s s n = s n+ s = s M . Next, we find that the set of integer numbers 
splits into four disjoint subsets: Z = A^i U N2 U N3 U N4 where Ni = 4Z + i. (We could 
have reduced solutions to only one subset by allowing < 5 < 4. However, higher order 
equations, such as order four here, naturally arise in loop quantum cosmology, and so we 
leave this explicit for the sake of generality.) 

When discretization is implied, the delta parameter is being fixed and the subset concept 
is very useful. If not, every subset gives rise to a channel Ni = N { + 5 with 5 G [0, 1). These 




(3) 
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four channels cover R completely. Equation (jHJ) can then be rewritten: 

4 +4 - (2 + 7 2 )^ + 4_ 4 = -A\n + 5\"s 5 n . 

For simplicity, we will drop the label 5 in s n and simply write s n . Different sequences will 
be distinguished from each other by the appearance of 5 inside the absolute value in 

Sn+4 - (2 + 7 2 K + Sn-4 = -A|n + S\ v s n . (4) 

3 Cosmological constant 

We start with the case of a cosmological constant, i.e. v — 1, in order to illustrate the main 
techniques. So equation (J3J) takes the form: 

s„+4 - (2 + 7 2 )s„ + s n _ 4 = -A|n + 5\s n (5) 

3.1 Generic asymptotic behavior 

We define 

S„ = (-1) 4 S n 

where n = 4m + q with g the subset index (q = i — 1 for TVj) and m an arbitrary integer. 
Introducing this into (0) results in 

s n+i + 2s n + s n -i = A\n + 5\s n (6) 

which is more convenient to investigate, starting with the n>l region. 
To establish the asymptotic behavior we define 

s n+4 =: h(n + A)s n . 

We have 

h(n + 4) = A\n + 5\-2-— - (7) 

h[n) 

which allows expressing h{n + 4) in terms of h(n). 

For n < 0, we can still keep n and 5 positive (n ^> A -1 ) but substitute n — > — n and 
<5 — > —5 in (jSJ). It is then more convenient to use s_( n+ 4)/s_ n instead, in order to investigate 
the regime of large negative n. This ratio is equal to 

S-( n +4) _ 1 



s- n h{—n) 

and the recurrence now is 

h(-n + 4) + , . l - = A\n + S\ - 2 . 
n{—n) 
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3.1.1 Continued fraction 

With ((Zj) and the analogous equation for h(n), h{n) = A|n + 5 — 4| — 2 — h{n — 4)~ 1 , we 
rewrite 

% + 4) = A\n + 5 1 - 2 j . 

A\n + 5-A\-2 



h{n - 4) 

This can be iterated, for which it is convenient to use a suitable notation to indicate that 
the k-th step of iteration has been performed: 

h (n + A) = A\ n + S\-2- 1 



h(n) 

/i x (n + 4) = A\n + 5\-2 



1 

An+5-4 -2 



/i(n-4) 

/i 2 (n + 4) = A|n + 5|-2 



A|n + 5-4| -2 

A\n + 5 



h(n 



h k {n + 4) = A\n + 5\ 



A\n + 5-A\ -2 



An + 5 



A|n + 5 - 4k\ 

By construction, 



h(n - 4A;) 



/i fe ( n + 4) = h(n + 4) for all fc = 0, 1, 2, 3, . . . . 

As we observe, in the k-th step we have k denominators with expressions as A| • | — 2 and 
an ending term —h(n — Ak)' 1 . This process of iteration allows investigating the behavior 
of h{n + 4) in terms of arbitrary initial conditions of the form s mo , s mo+ 4. 
For n < we perform the first step of iteration by 

1 1 

A|n + 5|-2 



h(-n) A\n + 5-4\-2-h(-n + 8) 

and have in the k-th step 

—j^- — - = A\n + 5\-2-— — (9) 

hW(-n) 1 1 A|n + 5-4| -2 v ; 

1 



A\n + 5 - Ak\ - 2 - h(-n + 8(fc + 1)) 
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with k = 1, 2, . 



3.1.2 Convergence 

The iteration process leads easily to asymptotic relations (large |n| behavior). In this 
section we will elucidate why this is the case, discussing only n > such that n = 4m, 
m 1. Then the ending term for the k-th iteration is —h(4(m — k))^ 1 . Also for simplicity 
we set k — m, and thus the two last denominators of h k {n + 4) are 



A|n + <S-4(fc-l)| -2- 



A|n + 5-4A;| - 2- 



h(n - 4k) 



(4 + 8)A - 2 + 



1 

2- A5 + 



HQ) 



To convince ourselves that choosing h(0) arbitrarily does not greatly affect the large n 
behavior of h(n + 4), we can simply reverse the whole procedure. We pick any value of 
h(0) and, after a few steps, the system looses any information about it. To see this we 
present some examples. Because we always have inverse powers of h(n), such as in 

h(4) = A5-2- 



h(0) 
and 

h(8) = (4 + S)A - 2 ' 



M4)' 

initial values for h(0) such that h(4) ~ (but not equal) could give rise to singular behavior. 
Indeed, h(8) is large, but then 

h(12) = (8 + 5)A - 2 - -L. ~ (8 + 5)A - 2 . 

So all information about the potentially dangerous value of h(0) has been lost and the 
iteration proceeds as in the generic case. 

For completeness, let us see now what happens when h(4) is equal to zero. To proceed 
we introduce a parameter e and a new initial value h e (0) such that 

lim h e (4) = AS - 2 - — ^- = and h £ (4) ^ for e ^ . 

As e tends to zero, we will find h(8) explode to infinity, and h(12) tend to 8A — 2. So this 
is just the limiting case of the previous one. 
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There exists, however, one very special choice of h(0) whose value is influential on the 
whole sequence h(n). If we pick h(0) such that the following sequence of an infinite number 
of constraints holds, 

- , 1 r ~ (4m + 5) A - 2 for all m G N , (10) 
h(4m) 

then cancellations occur within the continued fraction step after step such that the previ- 
ously observed "forgetfulness" cannot be realized. 

First of all, we can see that if the above sequence of constraints holds, it will hold only 
for one, unique initial value which we call /iboundcd(O) (unique for given A, S and k). In 
App. |X]we present more discussion on this choice and prove that it leads to a vanishing 
solution for large n. The discussion above and in App. [X] signals that bounded solutions 
are permitted by our asymptotic relations. Furthermore, we will present tools, which allow 
us to find bounded solutions in a simple and explicit manner. 

For values of h(0) slightly different from /ibounded(O), the iteration undergoes locally 
abrupt behavior, but after a few steps it smoothes again. So even for h(0) very close to 
^bounded (0) the information about the initial value is being lost. 

If h(0) is different from h hounded (0) it is useful to split the fraction (jSJ) into two pieces 
(recall that k — m — j). The lower part of the fraction, which is affected by initial value 
choices for which possibly the A| • | term does not dominate, will be called "part B". The 
upper part, unaffected by the initial value choice and for which A| • | — 2 is much bigger 
than one, will be called "part A" . Thus, there is a k-Q < k such that part B is identical to 
hk B (n + 4), and part A is hk-k B (n + 4) with part B used as initial condition. 

From our arguments, we conclude that part B will become a small initial value con- 
tributing to part A. So what is left is part A, which is convergent but gives a negligible 
contribution, compared to the leading A|n|, as we move to the large n regime. The same 
reasoning applies for any choice of an initial value h(n — 4k). It is also easy to convince 
oneself that the above arguments are independent of our previous choice n = Am. 



3.1.3 Generic growth 

If the initial value is different from ^bounded, there are at most finitely many near cancella- 
tions, and the following asymptotic relations (for large n) hold 

Sn+A 



s -(n+4) _ 1 



-h(n + 4) ~ 2 - A|n| 
2 - Aid 



s- n h(-n) 
since the fraction part of (JHJ or Q can be ignored. 



3.2 Pre-classicality 

The previous results allow us to draw conclusions for the pre-classicality of solutions, in 
particular alternating behavior. 
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3.2.1 Subset splitting 



The subset splitting presented in Sec. ED allows us to reduce 




— Alnls. 



n 



to four equations 



S4(m+1) 



2S4 m + S4( m _ i) 



A|4m + i + <5|s 4 , 



.to 



each valid on one set Ni (i = 1,2,3,4). Each of them is of second order, and thus needs 
two initial values. There are no mixing terms, i.e. terms that mix s n from different N 
subsets. So these four equations are fully independent. 

Because for large n (and for ratios different from Abounded) 



and analogously for s„(„ +4 )/s_ n , we see that the asymptotic behavior excludes solutions 
where s n and s n+ 4 (also s_ n and s_( n+ 4)) are of the same sign. This means that pre- 
classicality within one subset Ni is not preserved. Additionally, it may happen that the 
solution between different subsets is not pre-classical. Because, as we have seen, the original 
equation reduces to four independent equations of second order, for each equation we need 
two initial values. It may then happen, for example, that for n G N\ s n oscillates with 
different phase than for n e jV 2 . 

Since we have just proved that generic solutions of (0) are not pre-classical within a 
subset, non-preclassicality between different subsets is only of secondary interest. However, 
we will meet equations which are pre-classical within a subset, and therefore it would be 
of interest to remove non-preclassical behavior between different subsets on the whole set 
Z = Nx U N 2 U N 3 U N 4 . 

A simple example of such an equation is ©. The same reasoning as above implies 
now that s n+ 4 and s n (also s_( n+4 ) and s_ n ) are of the same sign, so the solutions are pre- 
classical within a subset. We will show how to remove non-preclassical behavior among 
different subsets by a proper choice of initial values. In the next part we will present 
generalization to other models. 

3.2.2 Determining regimes leading to solutions of constant sign 

Again, we start with the positive n direction of iteration which will simply generalize to 
the negative one. We rewrite © in a form 



s n+4 



2 - A\n 



s 



n 





where 



a(n) = A\n + S\ - 2 . 
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Shifting the label by four and substituting s n in the right hand side of (fTTj). we get s n+ 4 
in terms of s n _4 and s n _g. After k iteration steps we reach 

= f ik) (n)s n ^ k - g W (n)s n „ 4{k+1) (12) 

where / and g are functions of n whose structure is of course fc-dependent. For step k + 1 
we have 

JR? = f ik+1 \n)s n ^ (k+1) - g {k+1 \n)s n „ 4(k+2) 

and from 

s n - 4k = a(n - A(k + l))s n _ 4 (A;+i) - s n ^ 4 (k+2) ■ 
Substituting this in (fT2"j) we obtain 

fin* = (<*(" - A(k + l))f ik \n) - g {k \n))7s n _ A{k+1) - f {k \n)s n . A(k+2) . 

Comparing with previous expressions for s^+l we obtain the recurrence relations 

}< k+l \n) = a(n-A(k+l))f (k \n)-f {k - l \n) (13) 
( n ) = f( k )( n ) (14) 

having already used relation (fl^|l in (fT3*|) . To close these recurrence relations, we notice 
that 

/(°)( n ) = a (n) and g {0 \n) = /^(n) = 1 . 

We will now focus on relation (fT3*J) . which has a symmetric form similar to This 
suggests to introduce 

and rewrite equation (fT^j) as 

p(k + l,n) = a{n-4{k + l))- (15) 

p(fc,n) 



If we choose to use as initial values s 4mo+q and S4(m -i)+g (<? is the subset constant for N t 
defined by q = % — 1, and mo indicates the initial values in a given subset), we may write 



s n +i = f {k) (n)s 4mo+q - f {k 1) (n)s 4(mo _i )+ y (16) 

where n — Ak = 4m + q has been substituted by k = \{n — Am — q). To assure that, for 
example, the solution is positive for large values of n, the condition s n+4 > easily implies 

sgn(/)(s 4mo +g - r(^)s 4(mo _i) +g ) > 

with 

r(k) = r(mo,q) = lim 



p{k,n) 
9 



and sgn(/) is the sign of the / coefficient which is discussed in App. IB. II 

To find r(m , q) it is sufficient to use equation (JT5|) . Arguments on convergence of ()15|) 

are given in App. IB. 21 So by a proper initial value choice we can make all four solutions on 

Z positive or negative. If p is not discrete, one proceeds as before and makes all solutions 

positive or negative for all 5 6 [0, 1). 

Equation (JT5|) is easy to handle, but nevertheless it is instructive to discover some 

properties of this equation for a few particular cases (for simplicity putting 5 = 0). 



1. Subset Ni, s and s_ 4 as initial values: 
Then k = f and shifting variables in (j!5|) to 

p(k,n) = a(n - 4k) (17) 

a(n - 4(k - 1)) - - r rr 

v v a(n-A(k-2)) 

we obtain 

p{n/4, n) = -2 



4A-2 



3A - 2 

1 



nA - 2 



As expected, this expression is for large n almost independent of n. The only parameter 
left is A, of which two values are particularly interesting: 
a) A = 0: Because 

1 

= -1 



1 



-2 - 

we get 

lim 



n->oo p(n/4, n) 

in agreement with the results in App. IC.3I 

b) A ^> 1: On the right hand side of Eq. (|T7|) the fraction can then be neglected, and 
we have ^ 

lim — — ; = — . 

n-^oo p(n/&, n) 2 

This subset and these initial values are, however, very special because for k = k = j the 
A-term in a(n — 4k) in (|17|) does not contribute. In other subsets, this does not happen. 
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2. Subset N2, s% and s_3 as initial values: 

This time, n — 4k = 1 such that k = and 



p(k,n) =A-2 (19) 

5A - 2 



9A - 2 

We observe, as anticipated, the occurrence of A at leading order. This follows from the 
subset splitting where for subset N% we have 

S4(m+1) + 2s 4m + S4(m-1) = A|4m|s4 m . 

If we put m = in the above equation we see that the recurrence starts with S4 = — 2so — s_4 
which is A-independent. For N2, however, we have 

S4(m+1)+1 + 2S4, Tt +l + S4( m _i) + i = A|4m + l|s4 m+ i 

and putting m = leads to S5 = (A — 2)si — s_ 3 . So A appears at the zeroth iteration 
step. Now it is essential that due to (|T9*|) the regime for constant sign is sensitive to values 
of A in contrast to the subset iVj. It is clear that this is the generic behavior, Ni being the 
only exception. 

For negative n we proceed similarly, starting with s^/ n+4 % = /^(n)s_ n+ 4fc— g^(n)s^ n+4: ^+i) 
such that 

/( fc + 1 )(n) = a(n-A(k + l))f {k \n)-p- l \n) 
9 {k+1 \n) = f k \n) 

allows us to write s_ (n+4) = / (fc " 1) ( ra )^4(m -i)+q - f {k Kn)s 4mo+q with k = |(n + 4m + q). 
3.3 Properties of initial values 

In this subsection we assume positive n; the generalization to negative n is immediate. 

3.3.1 Asymptotic roles 

Using p we write 

S n+4 = f (k \n) (^Simo+q ~ ^7^T^(mo ( 20 ) 

and, using the ra-independence of p for large n, clearly see the different roles played by the 
two initial values. The leading order contribution to p for an arbitrary subset and arbitrary 
value of m is 

p(k, n) ~ A|4m + q + 5\ - 2 

and thus it can happen that p{k,n)~ l ~ is a small number, depending on where initial 
conditions are stated (mo) in relation to 1/A. From Eq. (J2U)) we conclude that the large n 
behavior of s n+4 then depends very weakly on s 4 ( mo _i) +(? . Thus, s imo+q dictates the rate 
of growth, which can easily be checked numerically as illustrated in Tab. ^ 



11 



S5 





1 


2 


10 


5 


5 


5 


5 




5 


5 


5 


5 





1 


2 


10 


10- I30 s 205 


147 


148 


148 


154 


3.42 


32.8 


62.2 


297 



Table 1: Dependence on initial values in the subset N 2 with A = 5, 5 = 0, mo = 2. 
3.3.2 Bounded solutions 

From the asymptotic relation it follows that the solutions are generically growing fast and 
in general are not bounded. However, using the fact (App. IB. 2)1 that for a(n) given in (jlljl 
the limit r(mo, q) exists we conclude that special solutions exist which are bounded. These 
are those solutions for which 

lim (s imo+q 7^— s 4(mo _i) +(? J = (21) 

i.e. 

Sim +q = f(mo, g)s4( mo -_i) +g . 

We easily see that r(m , q) = h hounded . This shows that the infinitely many conditions (fTUJl 
can indeed be satisfied and that the value Abounded exists. Most importantly, we now have 
a constructive procedure to compute the value for any given model via r(mo, q). 

If one puts s 4 ( mo _i) +g = and compares this with (jTUj) and the asymptotic relations 
one concludes that f( k W)(n) tends to (plus or minus) infinity and so is responsible for 
the asymptotic behavior. Putting s imo+q = one finds that the same must hold for 
/(*(»)-!) (n). Looking once again at (|16|) allows to understand why the solutions are 
bounded for Abounded and. why this value is unique (for given A, 5 and m ). Only then 
is the growth of S4mo+g^ (ji) canceled by the growth of S4(mo— i)+<j 

4 Arbitrary matter 

In this part, we investigate a general class of equations (@J): 

s n+A - b(n)s n + s n _ 4 = -A\n + 5\ v s n (22) 
with b(n) = 2 + 7 2 , and real v. We must distinguish two cases. 

4.1 v > 

Here, the generalization from the discussion of a cosmological constant is straightforward. 
First, we introduce 

s n = (-1) 4 Sn 
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(where q is the subset index). Equation (|2~2~j) then has the form 

s n +A + b{n)s n + s n _ 4 = A\n + 5\ v s n (23) 

We start with the positive n direction of iteration and, to obtain the asymptotic be- 
havior, define again 

h{n + 4) = — — 

which, from (|23jl . satisfies 

h(n + 4)=A\n + 5\ v -2-^--^-. (24) 

n[n) 

For the negative n direction of iteration we have 

S-( n +4) __ 1 

s_ n h(-n) 

and 

h{-n + 4) + — ^— = A\n + 5\ v - 2 - 7 2 . 

The convergence arguments of the preceding section still apply, and we obtain 

Proposition 1 Let s n be a solution to the difference equation A2ty) with v > 0. If the ratio 
of initial values is different from Abounded , then the rate of growth is given by 

Sn+A = _^+4 = + A)=2 + 1 2_ A | n |« (25) 

for large positive n, and by 

£_( !± 4) = _£ : > ± 4) = 2 + ^ _ A|n|B (2g) 

for large negative n. 

So we conclude again that solutions within a subset are not pre-classical. Possible 
non-preclassicality among different subsets can be removed applying a generalized regime 
finding method, presented in Sec. 14. 4[ to each subset equation (|23|). 

4.2 v < 

This model is very interesting physically, because it describes a universe with matter or 
radiation. As we will see, it contains many new features. Let us rewrite equation (|22jl in 
the form 

s n+ 4 - (2 + 7 2 - A\n + 5\ v )s n + s n _ 4 = 

and define 

F(n) = 2 + 7 2 - A\n + 5\ v . 

First, for negative v we have lim n ^ ±00 F{n) = 2 + 7 2 so we conclude that the results of 
Appendix will be useful, and Eq. (|4*5j) there is indeed an asymptotic equation. 



13 



4.2.1 7^0 

We use h(n + 4) = s n+ ^/ s n and reach 
h(n + 4) = 2 + 7 2 - A\n + 5\ v 



2 + 7 2 - A\n + 5-A\ 



2 + ^ 2 -A\n + 5-8\ v -- 



and 



j 2 + 7 2 - A\n + 5\ v - 1 



h{-n) 1 1 

2 + 7 2 -A|n + 5-4|^- 



2 + 7 2 - A\n + 5-8\ v 



In the large n limit, using results of IC.1.21 we have 

Proposition 2 Let s n be a solution to the difference equation Wty with v < and 7 7^ 0. 
// i/ie raizo 0/ initial values is different from h hounded , then the rate of growth is given by 

Sn+A ~-h{n + A) = 2 + 1 2 -x (27) 



s, 



for large positive n, and by 



for large negative n, where 



£ ^ = T7 1 ^ = 2 + 7 2 -x (28) 
s_ n h{-n) 



2 + 7 2 - V7 4 + 47 2 , 

X = r < 1 • 



We conclude that solutions are preclassical within a subset, i.e. s n+4 and s n (also s_( n+4 ) 
and s_ n ) are of the same sign for all large n (this means that also long distance oscillations 
are not possible). The only problem is to remove possible oscillating behavior between 
subsets, which is an easy task due to Sec. 14.41 (the limit of p(fc,n)~ 1 will exist) and is 
achieved by picking up proper initial values. 

However this does not exhaust all properties of the 7^0 case. Long distance and 
short distance oscillations are excluded for large n, but they still may occur for small n 
(the standard of smallness is parameter dependent). 

If one picks A and 7 in such a way that F{n) > 2 for all n, then (with help of Appendix 
EJ) solutions are always non-oscillating. But for other choices one in general has F(n) < 2 
for n e [—no, uq] and F(n) > 2 for n G Z \ [—no, no]. As long as F(n) < 2 the oscillating 
nature of solutions may manifest itself which, however, is not guaranteed: The condition 
F{n) = 2 has solution \ii\ = (7 2 /A) 1/Au = (A/7 2 ) 1 /^'. One may call n the value(s) of last 
oscillation. Depending on the phase of oscillation in a neighborhood of n, the solution 
may lose its oscillating nature just before or just after n. It is instructive to illustrate 
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-100 




Figure 1: Typical behavior of a solution with last oscillation at n = 1000 and small ttlq. 

this feature with a numerical example: A3, A = 5, 7 = 0.1, p = —1, S2 = 2, sq = 3. 
The solution passes zero for the last time at n = 474 while the value of last oscillation is 
n = 500 (see also Fig. |TJ). 

The necessary condition for small n oscillations is now easy to state: The value of last 
oscillation must be greater than 4mo + q or smaller than 4(mo — 1) + q, at least by a value 
of four. Otherwise, oscillations will not even have space to start. A sufficient condition 
will be formulated in what follows. 



4.2.2 7 = 

This case is profoundly different because F(n) < 2 for all n. Thus, with results in Appendix 
O we conclude that generic solutions are oscillating. We start investigations with the 
v — — 1 case, moving to the large n regime. For any e > we may always find such A<o 
that for A^ > Ao we have 2 — e < F{n) < 2. We can use (joTIj) in Appendix as an 
approximate solution, with one caveat. Since we are in the large n regime, we cannot pick 
arbitrary in ()5()j) . Rather, it should be some value at which the large n regime has 
already started. This means that, in general, we do not know anything about s no , s no _4, 
which are responsible for the amplitude of oscillation. However, the period of oscillation 
is independent of s no , s„ _4, and we may use (joT?|) to extract the period of oscillation. 
Using (}4T?|) and (joT?|) we find the oscillating part of the solution for large n to be given 
by sin(a;n/4) with period of oscillation T = 8tt/uj. It is more convenient to use half of 
the period, T := T/2 = Att/uj which will be the distance between two subsequent zeros of 
oscillation. (Because n is discrete, and its progress is four, it rarely happens that a solution 
takes zero value. Usually we have s n and s n+ 4 of different sign. By zero of oscillation we 
mean then a number z = n or z = n + 4, depending on which absolute value is smaller, 
s n or s n+ 4.) With u>(a) = tan~ 1 (v / 4 — a 2 /a) and substituting a by 2 — A/|n| we reach an 
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A 


initial values 


Z\ 


£2 


Z2 ~ Z\ 


f(n) 


1 


So = 2 S4 = 3 


332 


604 


272 


271.83 


2 


Si = 2 S5 = 3 


973 


1269 


296 


297.48 


10 


S10 = 2 S14 = 5 


1806 


1978 


172 


172.81 



Table 2: Numerical compared to analytical values for the half-period (1291) . 



approximate half period of oscillation of 



f(n) = r ~— -. (20) 

tan" 1 ((V4A|n| - A 2 /|n| 2 )/(2 - A/\n\) 



With tan 1 (z) = x — y + (x 5 ) for very large n we can write 

47T / 

T(n ~ -jfVWl (30) 
VA 

Since fl2llj) is indeed independent of s no and s no _ 4 , it should be a good approximation. 
To compare it with a numerical iteration we proceed as follows. Let Z\ and z<i be two 
subsequent zeros of a solution. We define (z± + z 2 )/2 = n and use this value in (J29*)l to 
check if T(n) ~ — 2^2 (• This is confronted with a few numerical results in Tab.0 showing 
indeed good agreement. 

As already mentioned it is not possible to derive an exact value for the amplitude 
from ([50)1 . However, it is possible to derive the rate of growth which is proportional to 
2(4A/|n| -A 2 /H 2 )~ 1/2 ~ 

Let us now drop the assumption v — — 1. For large n the oscillating part of the solution 
is sin(o;(n)n) with 



V 4A /i™ 


r - A 2 /|n 


2r 


2 - 


A/|n 


r 



^ 1 -1 

uj[n) = - tan I 

and r = |t>|. One concludes that oscillations may be suppressed when lim n _ +00 u(n)n = 
lim n ^ 00 exists which implies r > 2. So for 7 = we have 

Proposition 3 For —2 < v < £/ie solutions are oscillating with increasing period 
For u < —2 £/ie phase is decreasing so rapidly that oscillations are suppressed. The regime 
finding method presented in Sec. \4-4\ allows then to remove non-preclassical behavior among 
subsets. 

The oscillating nature of solutions for —2 < v < does not imply that they are not 
preclassical. The period of oscillation increases, and oscillations are of long distance type. 
For v close to —2, no oscillations occur for extremely large n. Fig. El illustrates the rate of 
growth of uj{n)n for v = —1.9. 

Using the results of this subsection, it is now easy to state the sufficient condition for 
small n oscillations in 14.2.11 The half period of oscillations for small n should be smaller 
than the value n — 4m — q for the positive n direction of iteration. The generalization to 
the negative n direction of iteration is again straightforward. 
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Figure 2: Phase u(n)n of the approximate solution with v = —1.9. 

4.3 Generic asymptotic growth 

Solving the difference equation for s n given by the asymptotic behavior we obtain a growth 
of the form A n//4 YYiLt K| u f° r v > from Prop. Qand c n with a constant c for v < and 
7 7^ from Prop. 121 For t> < and 7 = solutions are oscillatory with growing amplitude 
at large \n\. We summarize this by 

Theorem 1 (Generic asymptotic growth.) Letip^ be a solution to the difference equa- 
tion (QJ) without matter field. Ifv>0, then increases stronger than exponentially for 
\n\ — > oo. If v < 0, then |^| increases exponentially unless 7 = m which case the 
increase in amplitude is given by a power law depending on v. 

4.4 Initial value representation 

So far, we discussed mainly asymptotic properties where the precise form of the coefficients 
in (Q) is irrelevant. However, when we now turn to the initial value representation of 
bounded solutions, we must take into account that these initial values may be stated at 
small n where the asymptotic form of the coefficients is not precise. In order to prove that 
the initial value representation is still valid we thus have to use the most general form of 
our difference equation 

a(n)s n+4 - b(n)s n + c(n)s n _ 4 = f(n)s n (31) 

where b(n), c(n), f(n) are arbitrary functions, and for a(n) we only assume a(n) 7^ for 
all n even though generalizing to the case where a(n) has a finite number of zeros is not 
difficult. 

Naturally, this equation again splits into four subsets N±, N2, N3, N4. We make an 
ansatz 

s { nU = f {k \n)s n - ik - g {k \n)s n ^ {k+l) (32) 
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such that s!£+P = f^ k+1 \n)s n -^k+i) — g^ k+1 \n) s n ^^k+2) ■ On the other hand, from (flTT]) 

F(n-4(fc + l)) c(n- 4(fc+l)) 

Sn - 4fc " a(n-4(A; + l)) S ™- 4(fc+1) ~ a{n - 4(Jfc + 1)) Sn ~' {k+2) ^ 

with = b(n) + /(n). Defining 

. . F(n-A(k + 1)) , , 

h <"- 4 <* + 1 » = a (n-4(fc + l)) (34) 

and 

c(n — 4(fc + 1)) , \ 

n - 4 fc + 1 = ~) 777—^ 35 

a(n — A[k + I)) 

and inserting in (|32p using above definitions, we reach 

= /i(n - 4(Jfe + l))f {k) (n) - g {k \n) 
g^){n)=j{n-A{k + l))f k \n) 

and thus 

f k+1 \n) = h(n - A{k + l))f (k \n) - j(n - Ak)f^ k - l \n) . (36) 

Introduction of 

p{Kn) = WW) (37) 

shows that the initial value representation is still valid, and has the form 

mi \ ( ^'( 4m ° + 9) A /oox 

Sn+4 — J W I S4m +g p~(k~n) S4 ( m ° _1 ) +9 y ^ ' 

for the positive direction of iteration, the negative one following analogously. This repre- 
sentation is particularly important when the limit lim^oo p(k(n), n)~ l exists. Violations 
of pre-classicality among different subsets may then be removed and bounded solutions 
exist. We have 

Theorem 2 (Initial value representation of bounded solutions.) Leta(n), b(n), c(n) 
and f(n) be functions such that lim n ^ 00 p(k(n), n) _1 =: r(mo,q) exists, where p(k,n) de- 
fined in depends on the functions through < T%j) and (f^| ), 13*5)) F(n) = b(n)+f(n). 

Then the difference equation \31)) . such as (Qp ; in any subset N q+ i has a solution 
bounded at positive n whose initial values specified at mo satisfy S4 mo+q = r(mo, q)j(Am + 

mo— !)+?• 
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5 Example: Euclidean cosmological constant model 



For a Euclidean cosmological constant model the equation 

s n+4 - 2s n + s n _ 4 = Ans n (39) 

has been studied [9J, which is here of interest because special solutions are known (note the 
absence of absolute values around n in this toy model, which leads to explicitly known so- 
lutions). This can then be compared with our methods which easily apply to this equation, 
in particular the regime finding method and the condition for bounded solutions. 
For positive n 6 N\, Eq. (p?9*|) has solution 

Sn = cpj*.! (±) + c^y- + i (±) ■ (40) 



For n — ► oo the Y Bessel function tends to minus infinity and the J Bessel function to 
zero. Generically, the function Y thus dominates, which agrees with our results for the 
generic asymptotic growth (Theorem P): If v = 1, we have s n+ 4 ~ Ans n from Prop. ^ 
such that s n ~ (4A) n / 4 (n/4 - 1)!. With Stirling's formula n\ = V2imn n e- n (l + 0{n)) 
we obtain s n ~ 2\phx(n — 4) _1 / 2 (A(n — 4)) n / 4 e~ n / 4+1 . This can easily be seen to agree 
with the asymptotic growth of the Bessel function Y v {z) ~ — y/2/ 7rv(ez / 2v) ~ v in the case 
z = 1/2A, v = n/4 + z which determines the solution to the difference equation. 
If we choose s 4mo and s 4 ( mo _i) as initial values then C^ 1 is given by 

where det(m , A) = J mo+ ^Y mQ _ 1+ i_ - Y mo+ i_J mQ _ l+ i_ and a x = J mo+ _i_/ J mo _ 1+ _i_, a = 

•54m / s 4(mo-i)- ® ne can show that det(m , A) = det(0, A) > and J mo _ 1+ j_ f^r) > for 
all positive A and m . So if for example s^ mo ^i) +q > 0, then for a < a\ tlie solution for 
large n is negative and for a > ai it is positive. For a = a 1; C^ 1 = and the solution is 
bounded. In our previous notation (Theorem |2J), this means a± = r(mo, q) for the bounded 
solution. 

Comparing the results, we obtain a continued fraction representation for the ratio of 
values of Bessel functions, 

Jmo ~ 1+ ^ ( ^ } - 2 + 4Am 1 



J m +^(2A) 2 + 4A(m + l) 



2 + 4A(m + 2) - 



which in a similar form (with plus instead of minus signs in front of the fractions, which 
lead to modified Bessel functions of the first kind) can also be found in 

Generalization to other subsets and negative n is easy to perform, and also confirms 
the presented method. 



19 



6 Conclusions 



The methods introduced here can be used to determine the generic growth of solutions 
to a difference equation of type (J2J). Moreover, even in generically unbounded cases there 
are analytical conditions for the presence and values of special initial conditions giving rise 
to bounded solutions. These conditions can be evaluated numerically in a constructive 
manner, without being sensitive to numerical instabilities which would hamper any trial 
and error approach. We have illustrated this in Sec. El in a model where explicit solutions 
to the difference equation are known. 

Bounded solutions are often preferred, even though a clear reason is not known unless a 
physical inner product is used. One may impose boundedness as a final condition on allowed 
solutions, but usually in quantum cosmology initial conditions are imposed or arise at small 
volume. This, in general, does not guarantee that a wave function evolves to a bounded 
solution. In fact, bounded solutions are very special if there are classically forbidden regions 
at arbitrarily large volume. In such a case, growing solutions generically dominate. This is 
familiar from quantum mechanics, and has been shown here to be realized similarly at the 
level of difference equations directly. We emphasize, however, that the situation is more 
complicated with matter degrees of freedom, where the matter Hamiltonian is an operator 
and methods similar to those used here would require continued fractions of operators. 
Similarly, one has to deal with operator valued continued fractions in anisotropic models 
where now a(n), b(n) and c(n) as used in Sec. 14.41 are difference operators on labels other 
than n [THl I24j . Higher order difference equations can also occur for which we can 

still introduce h{n) as before. However, if there are also terms with s n ±%, the equation for 
h(n) obtained after dividing by s n contains s n+ s/s n = h(n + 8)h(n + 4) which makes it 
more complicated to solve. 

There are two different origins of unboundedness of solutions of difference equations. 
The first appears in classically forbidden regions, and is not different from unboundedness 
in the continuum limit. This happens if h(n) > 1, which can come from positive curvature 
7 > or negative matter terms such as a negative cosmological constant. There is also a 
second possibility, where now the behavior of solutions to a difference equation is completely 
different from that of a differential equation. This happens if h(n) < — 1, which we have for 
a positive cosmological constant. In this case, the solution is not only unbounded but also 
has alternating behavior and is thus not pre-classical. Such a behavior cannot occur in the 
continuum limit, and expresses the presence of large curvature. It should be interpreted as 
a breakdown of the minisuperspace approximation where the integrated curvature is the 
only curvature quantity. Local curvature is small if the cosmological constant is small, but 
the integrated curvature A|/x| becomes large just from large volume. As the analysis here 
shows, this behavior is clearly distinguished and easily detectable. 
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Appendix 



A The choice of h honnded 

Let us assume that the sequence of constraints (|10|) holds. This means that if we take a 
large positive value of n, then there exists a small ei(n) such that 

1 

h(n + 4) = A n - 2 - — — = eAn) 

h{n) 

If we assume h(n + 4) > h(n) then —h(n + 4) _1 > —h(n)~ l and 



+ 8) = A|n + 4| - 2 - — > A|n| - 2 - — — + 4A = 4A + e x {n) 

h(n + A) ~ V h(n) ' 



implying in turn 



h(n + 12) = A|n + 8| - 2 - — > A\n + 8| - 2 - 



1 



ft,(n + 8) ~ 4A + e^n) ' 

A cancellation is then blocked by the 4A term, which is in contradiction to the assumption 
of cancellations for all n. So we have h(n + 4) < h(n) and the sequence is decreasing. One 
can easily find out that h(n) with the choice Abounded is much smaller than one but positive. 
Together with the previous prove this implies that solutions are vanishing asymptotically. 

One can convince oneself, that this is the case by noting that the previous observations 
imply for large n (and positive solutions) < s n < c^ n-m -'' 4 s m for some m and < c < 1, 
and thus lim^oo s n = 0. For negative solutions we have > s n > o- n ~ m ^ A s m and again 
conclude lim^oo s n = 0. 

The generalization to large negative n is straightforward. 

B Sign of / coefficients and the existence of regimes 
B.l Sign of / coefficients 

By definition, the coefficients / obey the recurrence relation 

f k+l \n) = a(n - 4(Jfe + l))f {k) (n) - /^(n) 

with f(°>(n) = a(n) and /^ -1 '(n) = 1. We will prove that for a{n) > 2 for all n all / 
coefficients are positive. We have f^'{n) = a(n — A)a{n) — 1 > a(n) = f°(n) noting that 

a(n — 4)a(n) — 1 > 2a(n) — 1 = a(ri) + (a(n) — 1) > a(n) 

because a(n) — 1 > 1 > as assumed. Let us assume now that f®(n) > /^'^(n) for 
/ = 0, 1, . . . , m. Then 

f( m+1 \n) = a(n - 4(m + l))f {m \n) - f {m ~ l) (n) > 2f {m) (n) - / (m_1) (n) > f {m) (n) 
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and by induction we finish the prove that f (k \n) > 

quence f^ k \n) is growing strictly monotonically and its initial values are positive, implying 
f^ k \n) > 0. 

We now turn to the case where 

a(n) < 2 for n G [nj.,/^] and a(n) > 2 for Z \ [711,712] (42) 

and consider the positive n direction of iteration first. If we take large n (fulfilling n ^> 112), 
then a(n — 4k) > 2 for k = 0, 1, 2, . . . , fc where we have f^ k \n) > jn fc-1 '(n) > for 
k = 0, 1, 2, . . . , k , and so p(k, n) > for k — 0, 1, 2, . . . , Afo . We must then investigate 
what happens with the sign of / for k = ko + 1, . . . , k (one should note that k — ko does 
not depend on n). With help of the already established formula 

p(k + l,n) = a(n-4(k + l))-—^— (43) 

p(k,n) 

and p(0,n) = a(n) we find p(k Q ,n). We then use this formula again to see whether or not 
p(k + l,n) = /^ 0+1 ^(n)//^°)(n) changes its sign. If not, f( ko+1 \n) is also positive. We 
then investigate p(k + 2,n), and so on. Every time p(k,n) < we have a change of sign, 
until we reach the desired value p(k,n). So one such iteration procedure allows to find 
both p(k,n) and sgn(/) = / (S) (n) /\f^ k \n)\. 

We will prove in the next subsection that for a(n) as in (}4*2*J) lim n _ >0O p(fc(n),n) exists, 
and that p(k(n),n) is n-independent for large n, and so must equal p(k (n),n). So it is 
sufficient to apply the above procedure for one value of large n. 



B.2 Existence of the regime constant r(mo,q) 

The problem of existence of r(rriQ,q) = lim^oo p(k(n), n) -1 is closely related to the exis- 
tence of continued fractions. We will show this for the positive n direction of iteration. 
From Eq. (}4*3*j) we have 

1 

p(k, n) = a(n — 4k) 



1 

a{n -4(Jfe- 1)) 



a(n - 4(Jfe - 2)) 

Substituting k(n) — |(n — 4m — q) we reach 

p(k(n),n) = a(4m + q) . (44) 

a(4(m ° + 1) + " ) - a (4( m o + 2) + ,)-... 

So if a{n) is of the form (|42)1. there exists l such that for all I > l we have a(4(m +/) +q) > 
2. With help of results of IC.21 we conclude that (J44j) is convergent, and for large n almost 
n- independent. 
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C The case v = as an asymptotic equation 



If v < 0, we are dealing with a difference equation whose coefficients are almost constant 
at large n (or exactly constant for v = 0) and which is of the form 

s n+4 - as n + s n _ 4 = (45) 

with a real parameter a. This equation can be analyzed by explicit solutions, but care 
must be taken in interpreting the results since the limiting value a = 2 for large n and 
v < is a critical value of the difference equation. 



C.l Solutions 

We start with subset Ni and introduce the generating function F(x) = J2m=m s 4m^ m - 
From equation (|4*H|) we get (x _1 — a + x)F(x) — S4 mo x m ° _1 + S4( mo _i)X m ° = and 

F X = 1 1 2 ' 

1 — ax + x z 

To find solutions we must expand this generating function as a power series in x. One 
easily sees that the numerator will only result in a variable shift, and the main task is to 
expand 

G(x) = = V t n x n . (46) 

y 1 l-ax + x 2 ^ y ' 

n=0 

We will use 1 — ax + x 2 = (x — x\)(x — X2) with roots 



2 

which are real if and only if \a\ > 2. 



a + -y/a 2 — 4 a — -y/a 2 — 4 
xi = and £2 = (47) 



C.l.l \a\ > 2 

For this case, we write 

I = -!-P U 

(x — Xi)(x — X2) Ii - I2 \S - ^1 X — X2 J 

and (x — xi) _1 = — xf x (f — x/xi)^ 1 = —x^ 1 J2T=o( x / Xl ) k - F rom dlOj) we reach 

t n = (x 2 - n - x - x7 n - 1 ) = {a 2 - 4)- 1 / 2 (x^- 1 - x7 n ~ l ) . (48) 

X\ — x 2 

For the large m behavior we notice that for a > 2 we have x± > 1 and X2 < 1. For 
large n the second term in bracket in (J48|) becomes negligible and the first one results in 
lim^oo^ = 00. 
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C.1.2 \a\ < 2 

We make an analytic continuation to the principal sheet in ()48)1 and prove that the t n 
remain real. We see that X\ = \{a + iy/4 — a 2 ) = x* 2 and write x\ = \r\e tu} , x<i = \r\e~ tuJ . 
Before going further we notice |r| = 1 so that the roots are pure phase. We rewrite equation 
PI) as 

2 

t n = ==sin(a;(n+l)). (49) 

V4 — a 2 

Thus, for \a\ < 2 solutions are oscillating with constant amplitude and a = ao = 2 is a 
critical point. 

Fromcu = tan _1 (\/4 — a 2 /a) we find the oscillation period T = 2tt/lu. Since lim a ^2- T(a) = 
oo the closer a is to two, the longer the period is, and solutions slowly acquire the non- 
oscillating nature of \a\ > 2 fm lC.3l we show that also the case a = a = 2 is non-oscillating). 

Making the already mentioned variable shift in F(x) we reach 

Sim ^4mo^m—mo ^4(mo — l)tm— (mo+1) 

In terms of n and for an arbitrary subset we have 

Sn = S no tn-q_ mo — S no -4^2i^£_( mo + l) • (50) 

We also conclude that if a(n) is such a function that for all e there is an mo for which 
2 — e < a(n) < 2 for all n > m such that for large n we have 

a(n + 4) ~ a{n) ~ a{n — 4) (51) 

and (jBTH) with a a{n) is an approximate solution to s n+4 — a{n)s n + s„_ 4 = 0. 

It should be noted that the bigger the value of n is, the better condition (|5T|) is fulfilled. 
However, it does not imply that u>(n), obtained after substituting a by a(n), is constant. 
Rather, it varies very slowly compared to the change of s n when n is changed by four. 
One also must remember that a = = 2 is a critical point of the asymptotic equation, so 
substituting a by two instead of observing condition (JoTj) will not lead to correct results. 

C.2 Continued fractions 

We will employ another tool used to investigate loop quantum cosmology equations. This 
is a continued fraction, which we define by: 

1 

= x 

1 

a 

a — • • ■ 

Noticing (a — x)~ x = x we have 

x 2 - ax + 1 = . (52) 
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This is the denominator of G(x) in (j4T)j) and so existence of the continued fraction is 
related to the problem considered before. For \a\ < 2, Eq. (|52|) has no solutions and so 
the continued fraction does not converge. To see from the previous context why this is the 
case we introduce h(n + 4) = s n+4 /,s n in (jj^J) such that h(n + 4) = a — h{n)~ l . Thus, 

lim h(n + 4) = a 

a 

a — • • • 

and one sees the close relationship between equation (|45|) and continued fractions. For 
\a\ < 2 solutions are oscillating, and so is h(n + 4), so the continued fraction does not exist. 
For \a\ > 2, on the other hand, h(n + 4) = a — x and the continued fraction is well-defined. 

In particular, if we assume a > 2 then x = xi from ()47)1 . For finite n and if sufficiently 
many steps are taken, 

1 

1 

a 

1 



is almost ho independent, and so very close to x. 



C.3 The case \a\ = 2 

We start with the solution to (|4*H|) for a = 2, s n+4 — 2s n + s n _ 4 = which is 

s n +4 = (m + 2)s n _ 4m - (m + l)s n _ 4 ( m+ i) . 

We pick m in a such way that we have s n+ 4 in terms of initial values at n — 4mo and 
n - 4(m + 1), say: s n+4 = (m + 2)s n _ 4mo - (m + l)s n _ 4 (mo+l)- 

To obtain s n+ § in terms of the initial values, we shift the label from n to n + 4 and 
choose m such that ra + 4 — 4m = n — 4m , i.e. m = m + l. Thus, s n+8 = (m + 3)s n _ 4mo — 
(m + 2)s n _ 4 ( mo +i) implying 

"^n— 4mo 4(mo+l) ■ 

Because on the right hand side we have initial values, the progress of the sequence is 
constant. In the special case s n _ 4mo = s n _ 4 ( mo+ i), the progress is zero and solutions are 
constant. So solutions in the case a = ao = 2 are in general linearly growing (except for one 
choice of initial values). This should be compared with the oscillating nature for \a\ < 2 
and exponential growth for \a\ > 2. As we have already remarked, a = ao is a critical point 
of the asymptotic equation, resembling critical points of differential equations. 

It is also interesting to discuss how the general procedure to find bounded solutions 
works here. For the subset N\ we have s n+4 = (n/4 + 2)so — (n/4 + l)s_ 4 from which the 
condition 

n/4 + 1 

s ° > /a i o s - 4 
n/4 + 2 
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for a positive solution follows, and thus we must have so > S-4 in the limit n — > oo. Using 
our previous approach with the function p and noting that 




we have r(mo,q) = lim^oo p a =2(k(n), n) 1 = 1, reproducing s > S-4- 

The case a = — 2 follows immediately after the substitution s n = (— l)( n-9 )/ 4 s, 
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